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(Dated: May 3, 2006) 

We present a search for supersymmetry in the i?-parity violating resonant production and decay 
of smuons and muon-sneutrinos in the channels fi — > Xi M — y X2,3,4 Mi an d — > Xi % f 1 - We 
analyzed 0.38 fb _1 of integrated luminosity collected between April 2002 and August 2004 with the 
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DO detector at the Fermilab Tevatron Collider. The observed number of events is in agreement with 
the standard model expectation, and we calculate 95% C.L. limits on the slepton production cross 
section times branching fraction to gaugino plus muon, as a function of slepton and gaugino masses. 
In the framework of minimal supergravity, we set limits on the coupling parameter A211, extending 
significantly previous results obtained in Run I of the Tevatron and at the CERN LEP collider. 

PACS numbers: 11.30.Pb, 04.65. +e, 12.60.Jv 



Supcrsymmctry (SUSY) predicts the existence of a 
new particle for every standard model (SM) particle, dif- 
fering by half a unit in spin. The quantum number R- 
parity Q, defined as R = (_i)3S+l+2S ; where B L and 

S are the baryon, lepton and spin quantum numbers, is 
+1 for SM and -1 for SUSY particles. Often imparity is 
assumed to be conserved, which leaves the lightest super- 
symmetric particle (LSP) stable. However, SUSY does 
not require i?-parity conservation. 

If i?-parity violation (R p ) is allowed, the following tri- 
linear and bilinear terms appear in the superpotential Q : 

Wft p = ^X ijk LfL^E k + Xi jk LfQ^D k 

+\x'l jk U i i DiDl + l x i L i H l (1) 

where L and Q are the lepton and quark SU(2) doublet 
supcrficlds and E, U, D denote the singlet fields. The in- 
dices have the following meaning: i, j, k = 1, 2, 3 = family 
index; a, /3 = 1,2 = weak isospin index; £,ifi,C = 1, 2, 3 
= color index. The coupling strengths are given by the 
Yukawa coupling constants A, A' and A". The last term, 
/ij_LjiZi, mixes the lepton and the Higgs superfields. The 
A and A' couplings give rise to final states with multiple 
leptons, which provide excellent signatures at the Teva- 
tron. 

In the following we assume that only A 211 can be 
non-zero. This implies (muon) lepton number viola- 
tion. The R p coupling constants arc already constrained 
by low-energyexperiments, in particular A 211 < 0.059 ■ 
/lOOGeV |3|. For the squark masses rriq kincmatically 
accessible at the Tevatron, this limit on A' 211 is signifi- 
cantly improved by the present analysis. 

The DO Collaboration searched for resonant slepton 
production in Run I Q. The HI experiment at DESY 
searched for resonant squark production in the frame- 
work of i?-parity violating supersymmetry and published 
limits on the couplings A' ljfc . The combined limits from 
the LEP collider at CERN [||, assuming i?-parity vio- 
lating decay via LQZJ-couplings, arc m(x\)> 39 GeV, 
m(xf )> 103 GeV, m(^)> 78 GeV and m(fl)> 90 GeV. 

At pp colliders, an initial qq pair can produce a sin- 
gle smuon or muon-sneutrino assuming a non-zero A' 211 
coupling. The s-channel production is dominant and de- 
pends on the value of this coupling. The contributions 
of the t and u channels are negligible compared to the 
resonant s channel Q. 

The slepton can then decay into a lepton and a gaug- 
ino without violating i?-parity. The A 211 coupling allows 



neutralino decays via virtual sparticles (such as muon- 
sncutrinos, smuons and squarks) into two 1 st generation 
quarks and one 2 nd generation lepton. The \i decay 
branching fractions as predicted by mSUGRA, with the 
ratio of the Higgs expectation values tan/? = 5, the sign 
of the Higgsino mass parameter \i < 0, and the common 
trilincar scalar coupling A$ = 0, are assumed, leading to 
BR(x5 — > ^qiq'i) « BR(x? — * v^qiqi). In this analysis, 
the value of A 211 is always larger than 0.03, therefore the 
corresponding decay length is negligible. The dominant 
slepton intermediate decays as well as the corresponding 
final states are indicated in Table [I] 



I decay channel 


Dominant final states 


M -> xV 


2 /i, 2 jets 


M -> X^m 


1 /i, ]£t, 4 jets 




1 fj,, fir, 2 jets 


^ X ± V 


2 /i, 4 jets 



TABLE I: Smuon and muon-sneutrino decay channels: the 
final states correspond to xi — ¥ Mi<fi and xt ~ * ll'Xi- 

Because of the challenging multi-jet QCD environment 
and the advantage of the ability to reconstruct the neu- 
tralino and smuon masses, at least two muons were re- 
quired in the final state. This leaves the three channels 

(i) A XiM. (ii) P- X2,3,4M, and (iii) -> Xi,2 A* 
which are analyzed independently. The analysis is insen- 
sitive to events where the Xi decays into v^qq 1 and where 
no second muon is created in the cascade. 

The data for this analysis were recorded by the DO 
detector between April 2002 and August 2004 at a center- 
of-mass energy of 1.96 TeV. The integrated luminosity 
corresponds to 380 ± 25 pb" 1 . 

The DO detector has a central tracking system 
consisting of a silicon microstrip tracker and a central 
fiber tracker, both located within a 2 T superconducting 
solenoidal magnet, with designs optimized for tracking 
and vertexing at pscudorapidities 1 77 1 <; 3 and \r}\ < 2.5, 
respectively. A liquid-argon and uranium calorimeter has 
a central section covering pseudorapidities \r]\ < 1.1, and 
two end calorimeters that extend coverage to |7y| 4.2. 
The muon system covering |?7| < 2.0 consists of a layer 
of tracking detectors and scintillation trigger counters in 
front of 1.8 T iron toroids, followed by two similar layers 
behind the toroids. The first level of the trigger (level 1) 
is based on fast information from the tracking, calorime- 
try, and muon systems. At the next trigger stage (level 
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Invariant di-muon mass [GeV] Invariant di-muon di-jet mass [GeV] 



FIG. 1: Invariant di-muon mass in the two-muon sample (a) and reconstructed 4-body mass of two muons and two jets (b). 
The cascade decays in channels (ii) and (iii) lead to less energy per particle, thus lower invariant masses. The signal expectation 
for the point with m ; - = 260 GeV and m^o = 100 GeV is scaled in plot (a) by a factor of 100 and in plot (b) by a factor of 5. 
The dominant SM background is Z/'y* — > fifj,; other SM backgrounds are Z/y* — > rr; WW, WZ, ZZ, it and T production. 
The total SM uncertainty is shown as a dashed black line. The data is in good agreement with the SM expectation. 



2), the rate is reduced further. These first two levels of 
triggering rely mainly on hardware and firmware. The 
final level of the trigger, level 3, with access to the full 
event information, uses software algorithms to reduce the 
rate to tape to 50 Hz. 

The signal was simulated with SUSYGEN |||. The 
leading-order SUSYGEN signal cross sections have been 
multiplied by higher order, slepton-mass dependent 
Q CD-correction factors [H} of size 1.4 — 1.5 calculated 
with the CTEQ6M parton distribution functions 
(PDFs). The influence of the PDF uncertainty on the 
cross section is 3% — 6%, estimated from the CTEQ6M er- 
ror functions. The influence of the renormalization scale 
and the factorization scale hf is less than 5% for all slep - 
ton masses below 500 GeV, if m(i)/2 < fi F < 2-m(f) [l2 |. 

The dominant background is inclusive production of 
Z H* ~* W- It was simulated with the pythia 
Monte Carlo (MC) generator and normalized using the 
predicted next-to-next- to leading order cross section [J^ , 
calculated with the CTEQ6 PDFs. All other SM pro- 
cesses contribute only slightly to the total background 
as seen in Fig. ^ These contributions were simulated 
using the pythia and ALPGEN generators and normal- 
ized using next-to leading order cross section predictions 
calculated using CTEQ6M PDFs. All MC events were 
passed through a detailed detector simulation based on 
GEANT [l^ , followed by the reconstruction program used 
for data. 

Events were collected with di-muon triggers requiring 
at least two muons at level 1. At level 3 at least one 
track or one muon with a varying transverse momen- 



tum pt threshold of typically 5 — 15 GeV was required. 
To account for the trigger effects, simulated events were 
weighted using efficiencies determined from the data. 

All events were required to contain two muons. One of 
the muons was required to have pt > 15 GeV, and the 
second muon was required to have px > 8 GeV. A central 
track match was required for both muons. The muons in 
the signal are expected to be isolated. We define muons 
as "loose" ("tight") isolated, if the sum of the pr of the 
tracks in a cone with radius R CO nc = y A0 2 + Ar/ 2 = 0.5, 
where 77 = — In tan is the pseudorapidity and 9 is 
the azimuthal angle, around the muon direction is less 
than 10 GeV (2.5 GeV), and the sum of the trans- 
verse energies of the calorimeter cells in a hollow cone 
(0.1 < R conc < 0.4)) is less than 10 GeV (2.5 GeV). Both 
selected muons were required to pass the tight isolation 
requirement. The invariant di-muon mass distribution of 
this di-muon sample is shown in Fig. At least two 
jets with transverse momentum pr > 15 GeV and recon- 
structed with a cone algorithm (i? C onc = 0.5) 0] were 
required. Only jets within \rj\ < 2.0 were used. The re- 
constructed slepton mass with two muons and two jets is 
shown in Fig. ^p. The event selection is summarized in 
Table |n| 

Background from multi-jet QCD events was extracted 
from data using loose muon isolation requirements. This 
QCD enriched data sample was scaled to match the data 
in a signal free region. At least one isolation criterion 
with respect to other energy depositions in the calorime- 
ter or to other tracks must not be tight for at least one 
muon to separate QCD and the data sample. 
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FIG. 2: 95% C.L. limit on slepton production cross section times branching fraction to gaugino plus muon for the channels 
(i) ft — > Xi A* ( a )i (ii) A ~~ * X2,3,4 M (b) and (iii) — > A* ( c ) as a function of slepton and gaugino masses. The darkest region 
corresponds to a cross section of less that 2 pb. Successively lighter regions have successively higher limits. 
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FIG. 3: 95% C.L. exclusion contour on Ajn couplings within 
the mSUGRA framework for tan (5 = 5 and \i < 0. The arrows 
indicate limits on the slepton mass I, for a given coupling A211 . 



Two-dimensional selection requirements in planes 
spanned by the reconstructed I and x candidate masses, 
the invariant di-muon and di-jet masses, and the sums 



Cut 


Data 


SM expectation 


Signal eff. 


2[i selection 

p T jeti > 15 GeV 

Pt jet 2 > 15 GeV 


23206 
3852 
475 


22700 ± 70 ± 2900 
3760 ± 40 ± 560 
430 ± 10 ± 80 


5.5% ± 0.7% 
4.8% ± 0.6% 
2.4% ± 0.3% 



TABLE II: Expected and observed events at different stages 
of the event selection. The signal efficiency is given for the 
point with m ; - = 260 GeV and m^o = 100 GeV with respect 
to the total slepton production. The first uncertainty on the 
SM expectation is statistical, the second is due to systematics. 



of muon momenta and jet momenta were used to sepa- 
rate the signal s from SM backgrounds b. The selection 
requirements were chosen so that the signal efficiency x 
signal purity oc ^^ +b of a specific cut, applied on a train- 
ing sample, was maximized. The selection requirements 
were optimized for each (slepton mass, gaugino mass) 
combination (117 in total). 

In the fl —> Xi A* analysis (i), the slepton mass was 
reconstructed with the two leading muons and the two 
jets. In the signal MC, the leading muon usually orig- 
inates from the slepton decay vertex. The ncutralino 
mass was therefore reconstructed with both jets and the 
ncxt-to-lcading muon. 

Hadronic decays of vector bosons from the gaugino cas- 
cade to X? can lead to additional jets in channels (ii) and 
(iii). A simple likelihood was calculated for each combi- 
nation to reconstruct a vector boson and the neutralino 
candidate mass. The slepton mass was reconstructed 
from all jets with Et > 15 GeV and the two leading 
muons. 

After the optimization, for the point with mr = 
260 GeV and m^o = 100 GeV, we find 14/28/8 events in 
the data while 11.9±2.1±^ / 25.4±3.2^ / 6.b±1.6± 2 {° 2 
events are expected from SM backgrounds for the three 
channels, respectively, with a typical signal efficiency of 
up to 2%. For all 117 mass combinations, the data is 
in agreement with the SM expectation throughout the 
entire event selection range. 

The systematic uncertainties from different sources 
were added in quadrature-. For the limit calculation, the 
total systematic uncertainties of the background and sig- 
nal samples were taken to be 100% correlated. A sum- 
mary of the uncertainties is given in Table ITU with their 
contributions to the two muon and two jet sample. 

In the absence of an excess in the data, we set cross 
section limits on resonant slepton production. To be as 
model independent as possible, we calculated 95% C.L. 
limits with respect to the slepton production cross sec- 
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Uncertainty 


Background 


Signal 


Jet energy scale 


13.7% 


2 - 26% 


Muon ID 


7.8% 


8 -14% 


T.iiminnQitv (Hopq nnt annlv to Of^Dl 


5.5% 


6.5% 


Trigger efficiency 


5.2% 


4-9% 


MC a, A-factor, PDF 


3.7% 


5% 


QCD background estimation 


3.1% 




MC statistics 


2.2% 


3 - 24% 



TABLE III: Effect of the systematic uncertainties in the two 
muon and two jet sample on background and signal cross sec- 
tions. The muon ID contribution comprises the uncertain- 
ties due to muon reconstruction, isolation, track finding and 
matching, and resolution for the two muons. The systematic 
uncertainties on the signal strongly depend on the neutralino 
mass, so a typical range is given. 



tion times branching fraction to gaugino plus muon us- 
ing the Ch s method |l7j . The limit is then given in 
the slepton-mass and gaugino-mass plane, as shown in 
Fig. [5] In addition, our results are shown in Fig. [3] as 
A 211 exclusion contours interpreted within the mSUGRA 
framework, with tan/3 = 5, n < 0, and Aq = 0. The 
slepton-mass and gaugino-mass pair define the universal 
scalar and fcrmion masses mo and mi/ 2 . All three chan- 
nels were combined to form one limit for qq — > I, with 

I = fL,Vfx,. 

A lower limit on the slepton mass for a given LQD- 
coupling A 211 can be extracted from Fig. EJ These limits 
do not depend on other masses. They are indicated by 
arrows and summarized in Table IIVI Similarly, the ex- 
clusion contour can be translated within mSUGRA into 
constraints on other masses and parameters. 



Excluded slepton mass range 


Couplin 


I strength 


m(l) < 210 GeV 


for A211 


> 0.04 


m(l) < 340 GeV 


for A211 


> 0.06 


mil) < 363 GeV 


for A 2 n 


> 0.10 



TABLE IV: Limits on the slepton mass I for a given LQD- 
coupling A 2 n and tan/3 = 5, fj, < from Fig- El 

In summary, we have searched for i?-parity violating 
supersymmetry via a non-zero LQD-couplmg A 211 in fi- 
nal states with at least two muons and two jets. No excess 
in comparison with SM expectation was found and we set 
model independent cross section limits, improved com- 
pared to DO Run I by one order of magnitude. The lim- 
its are interpreted within the mSUGRA framework and 
translated into the best constraints to date on the cou- 
pling strength A 211 . DO Run I excluded slepton masses 
up to 280 GeV for X' 2U = 0.09 and m(xT) = 200 GeV. 
Now, slepton masses up to 358 GeV can be excluded, for 



A 211 ~ 0.09 independent of other masses. 
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